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ABSTRACT 

The timing capabilities and sensitivity of Kepler, NASA's observatory to find Earth-sized planets 
within the habitable zone of stars, are well matched to the timescales and amplitudes of accretion 
disk variability in cataclysmic variables. This instrumental combination provides an unprecedented 
opportunity to test and refine stellar accretion paradigms with high-precision, uniform data, containing 
none of the diurnal or season gaps that limit ground-based observations. We present a 3-nionth, 1 
minute cadence Kepler light curve of V344 Lyr, a faint, little-studied dwarf nova within the Kepler 
field. The light curve samples V344 Lyr during five full normal outbursts and one superout burst. 
Surprisingly, the superhumps found during superoutburst continue to be detected during the following 
quiescent state and normal outburst. The fractional excess of super hump period over the presumed 
orbital period suggests a relatively high binary mass ratio in a system where the radius of the accretion 
disk must vary by less than 2% in order to maintain tidal precession throughout the extended episode 
of superhumping. Disk radius is less restricted if the quiescent signal identified tentatively as the 
orbital period is a negative superhump, generated by a retrograde-precessing accretion disk, tilted 
with respect to the binary orbital plane. 

Subject headings: stars: novae, cataclysmic variables - stars: dwarf novae - stars: individual: V344 
Lyr - stars: white dwarfs 



1. INTRODUCTION 

NASA's Discovery mission Kepler (Borucki et al. 2010; 
Haas et al. 2010) was launched on 2009 March 6 upon 
a focused mission to detect and characterize terrestrial 
planets within the habitable zone of stars. The /sTep/er in- 
strument is a 0.95-m Schmidt camera with a 116 square 
degree field of view. The single photometric bandpass 
is broad, 95% of the response resides between AA 4,230- 
8,970A (Koch et al. 2010). Instrument commissioning 
performance is described by Caldwell et al. (2010). Ad- 
ditional to planet hunting, the photometric accuracy, ca- 
dence and uninterrupted integration of the same field 
over a multi-year mission provides unique opportunities 
to monitor sources of astrophysical interest. 

Keplers field of view contains ten known cataclysmic 
variables (CVs), listed in Table [1] with celestial coor- 
dinates, magnitude range, orbital period and sub-class, 
where known. These sources represent a specific and dy- 
namic stage in binary star evolution. The primary (more 
massive) stellar component is a white dwarf. In the past 
its progenitor star evolved off the main sequence and its 
expanding outer envelope engulfed a companion star. In 
a short-lived common-envelope episode, angular momen- 
tum was removed efficiently from the binary, leaving a 
close pair with separation a few stellar radii and orbital 
period a few hours (see Warner 1995 for a review). Sta- 
ble angular momentum loss continues slowly after the 
common envelope phase through magnetic braking and 
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gravitational radiation (Paczyhski 1967; Rappaport, Ver- 
bunt & Joss 1983). Through these angular momentum 
transfer mechanisms, the companion star's Roche lobe 
descends upon and makes contact with the companion. 
Contact initiates ballistic mass transfer from the com- 
panion star to the white dwarf through the inner La- 
grangian point, at which time the object is classed a 
CV. Mass transfer is stable and continues generally un- 
til the companion becomes degenerate. Thus, CVs have 
broad importance for constraining models of stellar evo- 
lution, binary evolution and the chemical enrichment of 
the Galaxy. 

The CV phase has another important feature. Ballis- 
tic material leaving the companion star cannot fall onto 
the white dwarf without losing angular momentum. In 
the absence of a powerful magnetic field from the white 
dwarf, material does this through the shearing layers of a 
circumstellar accretion disk (Shakura & Sunyaev 1973). 
Such disks occur in nature over a range of scales, from 
planet-building disks around protostars to those that 
feed the distant quasars. However, CV disks are unique 
in that their optical light dominates over both stellar 
companions, and systems are typically nearby, within 
a kpc. The mass transfer rate through the disk is of- 
ten modulated by a viscously driven limit cycle (Smak 
1984), resulting in intense, short-lived optical outbursts. 
Because of the observational detail available, CVs are 
primary targets for monitoring accretion disk behavior 
and testing accretion disk theory. 

Using the unprecedented detail of Kepler quiescent-to- 
outburst light curves we can test viscous accretion disk 
limit cycle models, quantify mass accretion rates, disk 
viscosity and trace the evolution of cooling and heating 
waves across the disk. Kepler will monitor the evolution 
of superhump periods throughout superoutbursts (Os- 
aki 1989), testing the viscous disk (Shakura & Sunyaev 
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Fig. 1. — The quarter 2 Kepler light curve of V344 Lyr. The insets are magnifications of the three boxes superimposed over the hght 
curve at discrete epochs. From left to right, insets contain a quiescent light curve dominated by photometric variability on, or close to, the 
putative orbital period, a quiescent light curve with no coherent orbital signal, and a superoutburst light curve dominated by double-peaked 
super humps. 



TABLE 1 

Identified cataclysmic variables in the Kepler, field. All data are extracted from the online catalog of Downes et al. 

(2001), EXCEPT where noted. 



Name 


KepID^ 


a (J2000) 


<5 (J2000) 


TTlv 




Sub-class'^ 


V344 Lyr 


7659570 


18 44 39.17 


+43 22 28.2 


13.8- 


-20.0 




SU UMa 


V358 Lyr 




18 59 32.95 


+42 24 12.2 


16.0- 


-20.0 




WZ Sge 


V447 Lyr 


8415928 


19 00 19.92 


+44 27 44.9 


17.2- 


-18.5 




dwarf nova 


MV Lyr 


8153411 


19 07 16.29 


+44 01 07.7 


12.1- 


-17.7 


3.18 


nova-like 


V452 Lyr 


7742289 


19 10 26.32 


+43:28:55.2 


17.6- 


-18.5 




dwarf nova 


V585 Lyr 




19 13 58.40 


+40 44 09.0 


14.9- 


-21.1 




SU UMa 


V516 Lyr 


2436450 


19 20 35.73 


+37 44 52.3 


18.9- 


-22.2 




dwarf nova 


V523 Lyr 




19 21 07.40 


+37:47:56.5 


17.7- 


-20.2 




VY Scl 




8751494 


19 24 10.82 


+44 59 34.9 


15 


.8 


2.94 


nova-like'^ 


V1504 Cyg 


7446357 


19 28 56.47 


+43 05 37.1 


13.5- 


-17.4 


1.67 


SU UMa" 



"Target identification number from the Kepler Input Catalog. 
'^Orbital period. 

"^See Warner (1995) for a definition of CV sub-classes. 
^Williams et al. (2010) 

1973), mass transfer instability (Osaki 1970), Enhanced 
Mass Transfer (EMT; Osaki 1985), Thermal-Tidal In- 
stability (TTI; Osaki 1989) and resonant disk (White- 
hurst 1988) models directly against unique photometric 
observations. We can also characterize disk flickering, 
quantifying its properties and the inferred mass transfer 
viscosity evolution as a function of luminosity and ac- 
cretion rate. Sensitive searches for the spin periods of 
the white dwarfs may reveal the magnetic nature and 
strength of the accretor and the truncation of the inner 
accretion disk due to magnetic pressure. We can also 
conduct systematic searches for orbital eclipses, yielding 
the inclination and mass ratio of the binary, greatly in- 



creasing the accuracy of accretion models. 

V344 Lyr is a dwarf nova which resides in the Kepler 
field. It is a faint source, with a quiescent magnitude of 
V ~ 19, and therefore has been observed sparsely until 
now. Optical photometry by Kato (1993) coincided with 
a bright epoch of V ~ 14 and revealed the presence of pu- 
tative 2.1948 ±0.0005 hour superhumps, indicating that 
this source is a member of the SU UMa class of dwarf 
novae. The normal outbursts occurred over a 16 ± 3 day 
cycle and the superoutburst recurred on a mean cycle of 
110 days (Kato, Poyner & Kinnunen 2002). Ak et al. 
(2008) used an orbital period-luminosity relationship to 
estimate a 619 pc distance to V344 Lyr, although the 
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Fig. 2. — Discrete Fourier Transform power spectra calculated over the duration of Q2 in 5 day slices, stepped by 0.5 days, 
minimum power, blue indicates maximum power. The putative orbital period, Porbi is detected at 11.7 cycles d~^ but is 
disappearing from the light curve before BJD 2,455,040. The putative supcrhump period, P^h, is detected at 10.8 cycles d 
confined to the superout burst, and is detected beyond BJD 2,455,080. 
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orbital period has, to date, not been confirmed directly. 

Relative to normal outbursts, superoutbursts are a 
magnitude or more brighter, considerably longer, have 
a shallower decay profile and reveal superhumps - pho- 
tometric oscillations on a period a few percent longer 
than the binary orbit (see Kato et al. 2009). The main 
competing models of superoutbursts, TTI and EMT, re- 
quire the accretion disk to expand in size up to or be- 
yond the radius where the rotation rate of the disk is 
in a 3:1 resonant state with the orbital period (White- 
hurst 1988). Additional fiux detected during superout- 
burst is the result of increased mass transfer through the 
accretion disk externally driven by the companion star 
through tidal dissipation. The superhumps are a combi- 
nation of tidally-driven precession of spiral arms in the 
outer disk and the reaction to precession of the classi- 
cal hot spot located where the ballistic accretion stream 
impacts the disk. 

During the first 90 days of observation, Kepler makes 
the unusual detection of superhumps within superout- 
burst, normal outburst and quiescent source states. This 
behavior is not predicted by the conceptual versions of 
superoutburst models. 

2. DATA ANALYSIS 

Kepler data are provided by the Science Operations 
Center to analysts in reduced and calibrated form after 
being propagated through the standard data reduction 
pipeline (Jenkins et al. 2010). Only the Simple Aper- 
ture Photometry (SAP) light curve provided for V344 
Lyr is considered in this paper. The Pre-search Data 
Conditioning (PDC) light curve, designed to optimize 
the search for 6 hour duration planet transits, is less well- 
suited to the accurate representation of high-amplitude 



photometric variability on timescales > 10 days. Con- 
sequently data artifacts that are removed in PDC data 
remain in our Kepler SAP. 

Artifacts include those caused by rapid attitude tweaks 
to correct pointing drift and those caused by the loss of 
fine-pointing. Fine-point loss is most often caused by 
thruster firings, performed on regular 3 day intervals to 
desaturate the spacecraft reaction wheels. Each event 
can affect Kepler data adversely for several minutes. Fur- 
ther gaps within the time series occur due to the space- 
craft entering occasional anomalous safe modes. Data 
are not taken during such events and recorded data fol- 
lowing a safe mode are often correlated with electronic 
warming over the first few days of operation. Further 
details of all artifacts and a list of their occurrence times 
are provided in the Kepler Data Release Notes #SQ. In 
this paper we make no attempt to correct data artifacts 
and simply remove them from the time series. 

The large panel in Figure [T] plots the short cadence (1 
min; Gilliland et al. 2010), barycentric-corrected light 
curve of V344 Lyi0 from the second quarter (Q2) of Ke- 
pler operation which spans the interval 2009 June 20 to 
September 16. Data gaps correspond to removed arti- 
facts and monthly data downloads. The quiescent V 
band magnitude is reported contemporaneously by the 
American Association of Variable Star Observers to be 
V ~ 19. While we make no attempt at a color correc- 
tion, these data indicate that Kepler short cadence data 
achieves 4% relative photometric accuracy upon V = 19 
sources, in agreement with Gilliland et al. (2010). The 
data reveal seven consecutive dwarf nova outbursts over 

^ http:/ /archive. stsci.edu/kepler/documents. html 

8 A Guest Observer Oflice target with KeplerlD: 7659570 
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the quarter of varying magnitude, duration and interval, 
with a comparatively short duty cycle on the order of 10 
days. One outburst is a superoutburst, confirming V334 
Lyr as an SU UMa class object, first reported by Kato 
(1993). 

Discrete Fourier Transforms over two short data sam- 
ples suggest we are detecting superhumps during super- 
outburst. Data extracted from BJD 2,455,006-2,455,025 
reveal a peak at what we will initially refer to as the or- 
bital period, Porb = 2.05737 ± 0.00001 hours, ahhough 
we argue for an alternative potential origin of this signal 
in Sec. [3l The error is a la estimate based upon 1,000 
Monte Carlo trials where the transform has been recalcu- 
lated after recasting each point in the photometric curve 
according to the inverse normal cumulative function and 
a random number generator. After 1,000 trials all pe- 
riod measurements occurred within the range 2.05733- 
2.05739 hours. A second epoch was extracted during 
the superoutburst over range BJD 2,455,058-2,455,080. 
While the power spectrum reveals a stronger but less 
coherent signal compared to the earlier epoch, the best 
period was determined by identical method to be Pgh = 
2.2023 ± 0.0001 hours and identified as a superhump. 
The fractional excess of the superhump period over the 
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orbital period is often characterized by e"^ 
-Forb)/^orb- If wc have identified the two signals accu- 
rately. The V344 Lyr e+ = 7.044 ± 0.005%. Conse- 
quently this binary would occur at the extreme end of 
the e+ distribution of the SU UMa class of stars tabu- 
lated by Patterson (2005). 

Neither the putative orbital nor superhump periods are 
permanent features of the Q2 light curve. We performed 
a dynamic Fourier power analysis of V344 Lyr by tak- 
ing 5 day segments stepped by 0.5 days and calculating 
Discrete Fourier transforms on each time slice. Power 
spectra are provided in Figure [2] The somewhat surpris- 
ing finding is that the quiescent photometric variability 
found in the post-superoutburst data at BJD 2,455,075 
does not have the same nature as the quiescent signal 
measured at BJD 2,455,010. Post-superoutburst, and 
throughout the following normal outburst, V344 Lyr con- 
tinues to be dominated by superhumps. The TTI model 
struggles generally to predict superhumps within quies- 
cence or normal outbursts and so the Kepler data provide 
us with an immediate challenge. 

3. DISCUSSION 

The detection of post-superoutburst superhumps in a 
CV is rare, but not unique. Two members of the SU 
UMa class, V1159 Ori and ER UMa have also revealed 
superhumps within quiescence or normal outburst (Pat- 
terson et al. 1995; Gao et al. 1999). The WZ Sge star 
EG Cue has also revealed superhumps up to 90 days after 
the decay from superoutburst (Patterson et al. 1998). 

Hellier (2001) has argued that quiescent and normal 
outburst superhumps are possible in extreme mass ratio 
objects where q < 0.07 (g — M2/M1; M2 is the compan- 
ion star mass. Mi is the white dwarf mass). Such a sce- 
nario is possible because angular momentum transfer by 
tidal dissipation is comparatively efficient in this regime; 
both cooling and heating fronts can propagate even when 
the outer edge of the accretion disk exceeds the 3:1 res- 
onance radius. Although no confirmed orbital period for 
this object exists currently, we ask in this section whether 
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Fig. 3. — The measured and Mi = 0.8 Mq zero-age main se- 
quence prediction of the orbital period-superhump period excess 
relation. V344 Lyr is represented by the red triangle. The inset 
plots the 3:1 resonance radius, Rs.i, and the tidal truncation radius 
of the disk, Disks precessing to the 3:1 resonance must reside 
between these two limits. Two vertical lines illustrate the mass 
ratio predicted by the model-dependent e'^-q relation of Patterson 
(1998) and the empirical relation of Patterson et al. (2005). 

V344 Lyr fits this same low-mass ratio profile. 

In order to proceed we must employ a small series of 
empirical relations. From Patterson (2005) we extract 
the orbital and superhump periods measured from a sam- 
ple of SU UMa stars. The distribution of e+ within this 
sample is plotted over Porb in Figure [3l Standard SU 
UMa objects are represented by circle symbols, while the 
three previous quiescent superhump examples are rep- 
resented by square symbols. Based upon the empirical 
fitting of superhump sources with relatively secure mea- 
sured mass ratios, Patterson (2005) provides the relation 
— q (0.18 + 0.299); we provide the predicted q on the 
right-hand vertical axis of Figure |3l The curve overlayed 
upon the distribution is a theoretical e~''-Porb prediction 
using the e~^-q relation and the assumptions that i) Mi 
= 0.8 Mq, and ii) the orbital period-companion star den- 
sity relation and zero-age main sequence (ZAMS) com- 
panion star mass-radius relation parameterized from the 
data of Clemens et al. (1998) are accurate; see Patterson 
(1998) for a more detailed description. The choice of Mi 
— 0.8 M0 is precipitated by a good, model-dependent 
match to the data in Figure [3] and the photometric tim- 
ing measurements of eclipsing dwarf novae which indicate 
a more massive white dwarf population in CVs relative 
to single white dwarfs (e.g. Littlefair et al. 2008). Fig- 
ure [3] reveals relatively low mass ratios for V1159 Ori, 
ER UMa and EG Cnc. 

The first of two scenarios we consider is that the high- 
frequency signal detected in Sec. 2 is the orbital period. 
In this case, the e+-predicted mass ratio of V344 Lyr, 
q = 0.27, is much larger than the other quiescent su- 
perhumpers. V344 Lyr is also not a good match to the 
Ml = 0.8 Mq relation; the ZAMS e+-g-Porb relation and 
it's dependent assumptions only remains good if Mi — 
0.58 M0. This is close to the mean mass of single DA 
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white dwarfs (Bragaglia, Renzini & Bergeron 1995), but 
would still clearly be a low-mass outlier within the mass 
distribution of CVs presented here. 

Within the upper- left inset of Figure |3] we plot the 3:1 
resonance radius, -Rsii, as a function of mass ratio and the 
tidal truncation radius, Rt, usually assumed to be 90% 
of the mean Roche lobe radius. Disks with outer radii 
residing between these limits will precess with the 3:1 res- 
onance mechanism. A disk in either the superout burst, 
normal outburst and quiescent states will potentially re- 
veal superhumps if it resides within this zone. Based 
upon the Patterson (2005) relation, this target should 
undergo precession but only if the outer radius of the 
disk remains within the range 0.444 < R/a < 0.454 dur- 
ing all three states, where a is the binary separation. This 
a relatively contrived parameter space and will provide 
an interesting challenge for the TTI and EMT models to 
reproduce the observed outburst phenomenology under 
such physical restrictions. Patterson (1998) provides an 
alternative e^-q relation, e+ — 0.23g/(l -I- 0.27q), based 
upon the argument that disk precession rate should be 
related to secondary star mass and that the most likely 
sites for tidal excitation are at disk radius R « 0.46a. 
This relation is of different form to that obtained empir- 
ically and provides a different mass ratio for V344 Lyr of 
q = 0.33. Under these circumstances, V344 Lyr falls even 
closer to the upper-g limit where superhumps are possi- 
ble, q = 0.34. This target could also undergo quiescent 
disk precession but only within a more highly contrived 
parameter space where the disk radius is very close to 
the tidal truncation radius throughout the superhump 
episode, 0.4368 < R/a < 0.4372. 

The second viable scenario we consider is that the 
shorter of the two signals detected is a negative super- 
hump (e.g. Wood, Thomas & Simpson 2009). This will 
occur if the accretion disk becomes tilted relative to the 
orbital plane, forcing retrograde tidally-driven preces- 
sion. Negative superhumps are comparatively rare. We 
note that of the 8 dwarf novae with putative negative su- 
perhumps compiled by Wood et al. (2009), the quiescent 
superhumpers VI 159 Ori and ER UMa are contained in 
that small sample. In such an instance, the binary or- 
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bit lies somewhere between our two measured periods. 
Making the assumption that Mi = O.8M0, then e+ lies 
in the range 4.6-5.3%. Consequently we can predict that 
the negative superhump fractional deficit lies within the 
range e" = 1.7-2.4% and that 2.086 < Porb < 2.101 
hours. The Patterson (2005) relation infers a mass ratio 
in this scenario of 0.20 < g < 0.22 and 0.453 < R/a < 
0.478 during superhump episodes. 

Both scenarios discussed can be tested with time- 
resolved optical spectroscopy that provides an indepen- 
dent measure of the orbital period. However both sce- 
narios infer that quiescent superhumps are not limited 
to extreme-low mass ratio CVs. 

4. CONCLUSION 

With its combination of photometric accuracy, time 
sampling and fixed field of view, Kepler opens a previ- 
ously unexplored window for the detailed monitoring of 
cataclysmic variables. Kepler reveals superhumps from 
V334 Lyr during quiescence, normal outburst and super- 
outburst. If the two detected signals within the light 
curve are the binary orbit and disk precession then the 
inferred mass ratio indicates that the outer disk radius 
varies by < 2% throughout multiple disk state changes, 
providing very restrictive physical limits for models of 
superoutburst to replicate. These restrictions become 
somewhat relaxed if we assume that the detected low- 
frequency signal is a negative superhump, in which case 
V344 Lyr is a target of high intrinsic value for accretion 
disk modeling. A radial velocity-derived orbital period 
will distinguish between these two scenarios. Ongoing 
monitoring with Kepler will also provide further indica- 
tions of the true nature and coherency of the observed 
signals. The continued, uninterrupted monitoring by Ke- 
pler oi this target, other known CVs, and newly identified 
objects in the field over the next 3 years or longer will 
provide a legacy for both the Kepler mission and CV 
community. 
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